dc-biased stationary transport in the absence of dissipation 
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We obtain stationary transport in a Hamiltonian system with ac driving in the presence of a dc 
bias. A particle in a periodic potential under the influence of a time-periodic field possesses a mixed 
phase space with regular and chaotic components. An additional external dc bias allows to separate 
effectively these structures. We show the existence of a stationary current which originates from the 
persisting invariant manifolds (regular islands, periodic orbits, and cantori). The transient dynamics 
of the accelerated chaotic domain separates fast chaotic motion from ballistic type trajectories which 
stick to the vicinity of the invariant submanifold. Experimental studies with cold atoms in laser- 
induced optical lattices are ideal candidates for the observation of these unexpected findings. 
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Bridging microscopies and macroscopics remains the 
ultimate goal in statistical mechanics, be it in or out of 
equilibrium. A considerable amount of work has been 
devoted to investigate this problem for low-dimensional 
stationary non-equilibrium states 0,0,1^. Stationary 
transport states under an external bias are often con- 
nected to nonzero dissipation (phase space volume con- 
traction and violation of time reversal invariance) , which, 
in turn, is equivalent to positive entropy production (for a 
recent review see Ref . jj] ) . Attempts to keep time rever- 
sal invariance by using Nose-Hoover thermostats 2] did 
not leave the grounds of dissipative systems and phase 
space volume contraction. 

The interface between microscopic dynamics and sta- 
tistical evolution and the role of nonlinearity and dissi- 
pation in nonequilibrium stationary states is an active 
area of research. Modern experimental studies provide 
an ideal testing ground to explore these problems using 
manipulations with cold atom ensembles in optical poten- 
tials 5]. The driven pendulum, a paradigmatic system 
for the study of dynamical chaos [Sl, is realized using a 
periodically modulated optical standing wave, with the 
possibility to control the strength of dissipation down 
to arbitrarily small values For example, recent cold 
atoms experiments 7] study the crossover from dissipa- 
tive to Hamiltonian ratchets, confirming theoretical pre- 
dictions 0. 

In this Letter we explore the route of obtaining station- 
ary transport using a dc bias in the absence of dissipa- 
tion and corresponding phase space volume contraction. 
A necessary condition for such a stationary, bias driven 
current to occur is a mixed phase space l6| with coex- 
isting regular invariant manifolds and chaotic regions for 
the unbiased case. A consequence of this coexistence is 
that directed transport may arise locally on regular in- 
variant components of phase space. While the chaotic 
phase space regions transform into accelerated evolution 
for any arbitrary small value of the dc bias, the regu- 



lar invariant submanifolds (unstable periodic orbits, reg- 
ular islands and cantori with nonzero mean velocities), 
initially embedded in the chaotic region, persist under 
a finite dc bias. Moreover unstable periodic orbits at 
zero dc bias lead to the appearance of stable islands for 
some finite nonzero value of the dc bias. The presence 
of cantori 9] and the corresponding sticking of chaotic 
trajectories in the vicinity of regular manifolds is man- 
ifested by the chaotic phase space part which, though 
accelerated, shows large time separation from chaotic tra- 
jectories which are initially located far from the regular 
structures. 

We consider the canonical Hamiltonian model of a 
particle moving in a one-dimensional space-periodic po- 
tential U{X) = ^ cos(27ra;) under the influence of a 
time-periodic space-homogenous external field E{t) = 
EacSin{ujt) 6], to which we add an external dc bias Ei,: 



X = p, p = sin(27ra::) -|- Eac sm{ujt) — Ei, 



(1) 



Due to time and space periodicity of the system we can 
map the original three-dimensional phase space {x,p,t) 
onto a two-dimensional cylinder, = (x mod\,p), by 
using the stroboscopic Poincare section after each period 
T = 2tt/uj, cf. Fig. 1(a). Although the dynamics is peri- 
odic in X direction, it contains the complete information 
about transport in the extended system (1). The velocity 
X — p along a trajectory is the same for both cases. 

In the case of zero dc bias, Et = 0, the phase space of 
the system is characterized by the presence of a stochas- 
tic layer which originates from the destroyed separatrix 
of the undriven case Eac = 0. The chaotic layer is 
confined by transporting KAM-tori, which originate from 
perturbed trajectories of particles with large kinetic ener- 
gies. These tori are noncontractible (since they cannot be 
continuously contracted to a point on T^), and separate 
the cylinder. 

The stochastic layer is not uniform and contains regu- 
lar invariant manifolds. There exists a whole hierarchy of 



embedded regular islands, which are filled by contractible 
tori, cf. Fig. 1(a). The centers of the islands contain el- 
liptic periodic orbits XTp(i) = {xTp{t),PTpit)} with the 
period Tp = kT, k = 1,2, ... and integer shifting distance 
L, such that 

XTAt + Tp)=XTp{t) + L, pTp{t + Tp)=pTp{t) (2) 

The transport properties of such an orbit (and of all tra- 
jectories inside the corresponding island) are character- 
ized by a winding number v = L/Tp. If w ^ then 
that PO and the corresponding island are transporting. 
In addition, within the chaotic sea there is an infinite 
number of unstable periodic orbits (UPO) with different 
periods and winding numbers. UPOs are unstable with 
respect to small perturbations and they are not isolated 
from the chaotic layer. Finally cantori exist which 
generate semipenetrable barriers in phase space. 
Eq.l^ is invariant under time reversal 



S -.t^ -t + T/2, 



X, p 



(3) 



which changes the sign of the current J — v. Thus, trans- 
porting invariant POs and islands appear as symmetry- 
related pairs in phase space with opposite winding num- 
bers. 

Contrary to the common expectation that all trajecto- 
ries acquire unbounded acceleration for nonzero dc bias 
-Ef, > 0, recent numerical evidence tells that this is not 
true for all trajectories ^^l- Fig- 1(b) we plot the 
escape time Tesc to be accelerated below the threshold 
Pt = —10 as a function of initial conditions {xo,Po} 
for nonzero E^. The window {xo,po} coincides with 
Fig. 1(a). The escape time is reached by a trajectory if 
p{Tes) = —10, and we integrate up to t = 200. While 
most of the chaotic layer in Fig. 1(a) is quickly acceler- 
ated, regular islands persist and their trajectories are in 
fact not accelerated at all. Moreover we observe phase 
space regions with delayed acceleration. These trajecto- 
ries stick for large time to the persisting regular islands. 

Noncontractible KAM-tori do not survive in the biased 
Hamiltonian system. If there exists at least one accelerat- 
ing trajectory, then KAM-tori do not persist Such 
accelerating trajectories always exist for Ef, ^ 0, both 
with Eac = and Eac 7^ 0. 

Let us give analytical proof that POs at Et — persist 
for nonzero E^. In order to formulate the proof we define 
the stability property of a PO (t). Towards this goal 
we linearize the phase space flow around it, and map it 
onto itself by integrating over one period Tp. The result- 
ing 2x2 symplectic Floquet matrix M has ei genv alues 
(Floquet multipliers) Ai and A2 with A1A2 = 1 [IjI- For 
a stable PO, both multipliers are located on the unit cir- 
cle, while for an UPO, both multipliers are located either 
both on the negative or positive real axis. 

Denote a solution of with initial conditions 
Xo = {xo,po} at t = by X{t,'Ko,Eb) = 
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FIG. 1: (a) Poincare section for the system Eq.(l), uj = 
2tv,T = l,Eac = 5.8, for (a) Et, = and (b) Et = 0.183. 
Several POs are shown: one stable {k — 1, L = l)(squares), 
and several unstable POs {k — 2, L — 1) (circles), (fc = 3, L = 
2) (upward triangles), and {k — 3,L = 1) (downward trian- 
gles). The island group near p — has zero winding number 
V = 0. (b) Time Tesc to be accelerated below the threshold 
Pt = — 10 as a function of the initial conditions in phase space. 
Trajectories from the islands are not accelerated at all. Inset 
in Fig. 1(b): Poincare section zoom of phase space structure 
showing the transformation of an UPO into a stable PO with 
a surrounding invariant regular island. 



{x{t,xo,Po,Eb),pit,xo,pQ,Eb)}. If Xo is close enough 
to Xtp (0) and Et is sufficiently small, then the solution 
'X.{t, Xq, Eb) on the interval [0,Tp] is close to the unper- 
turbed PO. X is Tp-periodic if and only if the vector 
function 



F(Xo,i?fc) 



x{Tp, xo,po, Eb) - xq - L 
p{Tp,xo,po,Eb) -po 



0. (4) 



Since X(t, (0),_pTp (0), 0) = Xxpit), Eq. ^ is satisfied 
by Xt (0) at Eb = 0. By the Implicit Function Theo- 
rem PJ, if the Jacobian of the vector function F(Xo, Eb) 
with respect to Xq is non-zero at Xt'p(O), Eb ~ 0, 
then the solution Xtp(O) can be continued to Eb ^ 0, 
Xo(i?6) {xQ(Eb\pvi{Eh)\ which corresponds to a Tp- 
periodic orbit of the perturbed system. The Jacobian of 
the function ^ at Xtp(O) is det9F/aX = det[M - I], 
where I is the unity matrix. The Jacobian is nonzero if 
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FIG. 2: Variation of Floquet multipliers, \j = |Aj|e'*J, of 
PO's vs bias Eb for Eq.(l) with parameters and symbols as 
in Fig.(l). We show both multipliers for stable PO (squares) 
and I All > 1 for unstable POs. 

^^1,2 7^ 1 for M. Thus, almost all POs of the unperturbed 
system persist for nonzero dc bias since generically their 
Floquet multipliers are not equal to unity. POs can be 
continued up to a value of Ei, where their Floquet mul- 
tipliers collide at +1. Linearly stable POs in Hamilto- 
nian systems are always enclosed by quasiperiodic tori, 
which form a regular island. Thus, the transporting con- 
tractible islands also persist for nonzero dc bias. Note 
that the symmetry Q persists for nonzero Eb. Conse- 
quently any invariant transporting manifolds which per- 
sist for nonzero Eh come in pairs and stationary currents 
occur in both directions 

We expect that for large enough dc bias i?b ^ 1, a 
transporting PO will disappear. This occurs when its 
Floquet multipliers collide at -1-1, X±(E^) = \2{El) = 1. 
This is a saddle-center bifurcation when two POs of 
the same period coalesce, one of them being initially sta- 
ble and the other one unstable. Since stable POs are en- 
closed by regular islands, this can happen only when the 
corresponding island surrounding the stable PO shrinks 
to zero. 

In Fig. 1(a) we depict the phase space structure of 
Eq.lUJ for w = 2Ti,Eac = 5.8 and Eb = 0. Using the 
Newton method, cf. e.g. 0|, we find stable and unsta- 
ble POs (various symbols in Fig.l), and continue them to 
nonzero values of E^- Stable POs exist up to Ei, « 0.412. 
The evolution of the corresponding Floquet multipliers 
A = |A| exp(i0) is shown in Fig. 2. 

Transporting UPOs become more stable with increas- 
ing dc bias, since they have to acquire Ai,2 = 1 at a crit- 
ical dc bias value before disappearing. Moreover, UPOs 
with negative real Floquet multipliers have to become 
stable POs before a critical bifurcation happens, since 
the only way to move the multipliers from the negative 
real axis to -1-1 is to transport them to -1 and subse- 
quently to shift them along the unit circle. During that 
last stage an originally unstable PO transforms into a sta- 
ble PO, and a new transporting regular island emerges 
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FIG. 3; Solid line: spatial distribution A'^(a;) after t = 100 for 
LJ = 27r, Eac ~ 5.8 and E-b — 0.183 and when 2.5 ■ 10* trajecto- 
ries started out from the uniform distribution in the rectangle 
(x £ [0.3, 0.75], p G [0.7,1.4]), which enclose the transport- 
ing island with v = \. The upper inset shows initial condi- 
tions which correspond to the tail of N{x), x £ [—200,50]. 
Lower inset depicts the Poincare section of a single trajec- 
tory that started out in the "sticky" vicinity of the island. 
Dashed line: distribution N(x) with initial conditions shifted 
by Apo = —2.2 and additional shift in x (see text). 

in phase space. Thus, along with an overall shrinking of 
the size of transporting regular islands with increasing 
Eb, new islands are born as well. In the inset in Figl(b) 
a zoom of the Poincare section shows such a case. For 
Eb = we found an UPO with Tp = 2 and winding num- 
ber V = 1/2 and Ai w -8.365, A2 = 1/Ai « -0.119. This 
orbit becomes stable at Eb ~ 0.178 (when the multipliers 
collide at —1 and enter the unit circle) and disappears at 
Eb « 0.185, when Ai = A2 = 1. In between these two 
bias values the now stable PO is surrounded by a regular 
island of finite size. 

For any PO the momentum p(t) is a periodic function 
and can be represented in the form p{t) = v +p'^(t), such 
that p"it + Tp) = p"{t), (p°(t))Tp = 0. Then we obtain 
the energy balance equation: 

vEb={p%t)Eac{t))T^. (5) 

An external ac-field mimics an additional effective fric- 
tion in the nonequilibrium Nose-Hoover oscillator |l5| : it 
pumps energy into the system in order to compensate the 
work against bias (for POs with positive v) and acts as 
an energy sink to compensate the acceleration in the case 
of current along a bias (for POs with negative v). But 
in our case, at variance with the Nose-Hoover system, 
the generating equations are Hamiltonian and dissipa- 
tionless. 
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All trajectories except the regular islands and the POs 
will be accelerated without bound for nonzero E},. Nev- 
ertheless, the mixed phase space structure for Et, — 
leads to a strongly nonuniform acceleration in time for 
trajectories with different initial conditions. That can be 
already observed in Fig. 1(b). Indeed, for Et ~ chaotic 
trajectories in the vicinity of regular islands may stick to 
the island boundary for long times and produce a unidi- 
rectional flight with the velocity v due to the stickiness 
effect ^3 . The origin of such a behavior is the presence 
of partial barriers in phase space formed by cantori . 
This feature is also observed on finite times for nonzero 
Eh. Trajectories may stick for long times to regular island 
boundaries without acceleration (see Poincare section in 
lower inset in Fig. 3). Only after some rather large escape 
times Tesc such a trajectory will eventually leave the is- 
land vicinity and suddenly accelerate. In Fig. 3 (solid 
line) we compute the distribution N{x) of displacements 
of trajectories with initial conditions uniformly covering 
a phase space part which includes a transporting island. 
We observe that the transient time T^sc can be extremely 
long for initial conditions in the island's vicinities (see 
also Fig. 1(b)). In Fig. 3 the right peak corresponds to 
initial conditions inside the island (no acceleration for all 
times). The left front peak corresponds to initial condi- 
tions that are accelerated right from the beginning. In 
the absence of trajectories which stick to regular islands 
with a broad distribution of escape times, the two peaks 
would be separated by a depleted region, i.e. a gap. In- 
stead we observe a smeared out distribution which corre- 
sponds to trajectories with long transient sticking to the 
island boundary. Indeed, initial conditions which corre- 
spond to the gap region stick to the island (upper inset 
in Fig. 3). A shift of the initial conditions by Apo = —2.2 
down into the homogenous part of the phase space in 
Fig. 1(b) yields a distribution N{x) which contains a sin- 
gle peak due to immediate acceleration (Fig. 3, dashed 
line). Note that we shifted N{x) here by Axq — 220 
which occurs for particles with V{x) = and Eac = if 
the initial momentum is shifted by Apo — —2.2. 

In conclusion, we show that stationary, dissipation- 
less transport is sustained by a nonzero dc bias. The 
dc bias provides a unique possibility to separate the dif- 
ferent parts of a mixed phase space, since it accelerates 
the chaotic part away and leaves the regular invariant 
manifolds basically untouched. Moreover, even within 
the accelerated part, proper distribution functions after 
finite acceleration times reveal the intriguing properties 
of surviving cantori and sticking. We give an analytical 
proof for the persistence of periodic orbits and thus for 
regular islands. We found evidence for the persistence of 
cantori as well. Our results are not limited to the case in 
Eq.(^ and are valid as well for the general case of a po- 
tential U{t, x) = U{t + T,x) = U{t, x + L). The presence 



of weak dissipation will not change the situation drasti- 
cally. Since it can be treated as a smooth perturbation, 
POs will persist for nonzero damping. Stable POs be- 
come attractors (limit cycles), and unstable POs saddles 
[l^ . A side result is then that an infinite number of sta- 
ble POs for zero damping will lead to an infinite number 
of attractors for infinitesimally weak damping |l7|. 

A possible experimental realization can be obtained 
with cold atoms in a standing wave potential , where an 
additional dc-bias can be realized by a gravitational force. 
In such a case, all of the above results could be verified. 
The persistence of transporting regular islands and POs 
can be used for the preparation of monochromatic matter 
waves. 
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